Initial state geometry and the role of hydrodynamics in proton-proton, 
proton-nucleus and deuteron-nucleus collisions 



o 



Adam Bzdak,^ Bjorn Schcnkc,^ Prithwish Tribedy,'' and Raju Vcnugopalan^ 

^RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 
^Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 
^ Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata 700064, India 

We apply the successful Monte Carlo Glauber and IP-Glasma initial state models of heavy ion 
collisions to the much smaller size systems produced in proton-proton, proton-nucleus and deuteron- 
nucleus collisions. We observe a significantly greater sensitivity of the initial state geometry to details 
of multi-particle production in these models compared to nucleus-nucleus collisions. In particular, 
we find that the size of the system produced in p-|-A collisions is very similar to the one produced 
in p-|-p collisions, and predict comparable Hanbury-Brown-Twiss radii in the absence of flow in 
both systems. Differences in the eccentricities computed in the models are large, while differences 
amongst the generated flow coefficients V2 and V3 are smaller. For a large number of participants in 
proton-lead collisions, the V2 generated in the IP-Glasma model is comparable to the value obtained 
in proton-proton collisions. Viscous corrections to flow are large over characteristic lifetimes in the 
smaller size systems. In contrast, viscous contributions are significantly diminished over the longer 
space-time evolution of a heavy ion collision. 
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INTRODUCTION 

A recent discovery in high multiplicity proton-proton 
and proton-nucleus collisions are correlations between 
pairs of charged hadrons that are collimated in their 
relative azimuthal angle and are long range in relative 
rapidity. These "ridge" correlations were mostly unan- 
ticipated for such small size systems, albeit a similar 
striking effect was previously seen in heavy ion colli- 
sions at RHIC and subsequently also at the LHC. The 
ridge in ^/s = 7 TcV proton-proton collisions was dis- 
covered by the CMS collaboration [l|. In proton- lead 
collisions at y/s = 5.02 TeV/nucleon, a sizable ridge was 
observed by the CMS coUaboration [H, the ALICE col- 
laboration [3| and the ATLAS collaboration [J|. In ad- 
dition, the PHENIX collaboration at RHIC recently dis- 
covered a ridge in very central deuteron-gold collisions at 
yi = 200 GeV/nucleon @. 

These long range correlations are of fundamental im- 
portance because they probe the very early time dynam- 
ics of matter produced in hadronic collisions. A question 
of considerable recent interest is whether the ridge effect 
in p+p and p/d-|-A collisions is due to initial state effects 
arising from the correlations of gluons already present in 
the nucleon and nuclear wavefunctions, or whether it is 
due to final state rescattering effects that are amenable to 
a hydrodynamic description. In both cases, one assumes 
a dynamical scenario where long range rapidity correla- 
tions arc generated in the initial state. [J The question is 



whether the azimuthal collimation observed in the ridge 
is also due to the same initial state correlations that gen- 
erate long range rapidity correlations, or whether they 
are generated primarily by the final state flow of these 
correlated structures. 

A powerful framework in which long range rapidity 
correlations can be computed systematically is the Color 
Glass Condensate (CGC) effective field theory (EFT) ^. 
In the CGC EFT, these correlations are a consequence 
of gluon saturation at central impact parameters in the 
proton and nuclear wavefunctions. In the hadronic col- 
lision, gluon fields are generated that stretch out in ra- 
pidity between the receding hadrons, and are coherent 
in the transverse plane over distances 1/Qs, where Qs is 
the saturation scale. The saturation scale in a hadron or 
nucleus is a function of the parton momentum fraction x 
and impact parameter, and grows with increasing energy 
and nuclear size. 

Multi-particle production, by the decay of the gauge 
field configurations corresponding to these Glasma flux 
tubes [7[ , is nearly boost invariant and nearly azimuthally 
isotropic; the resulting multiplicity distribution is the 
negative binomial distribution [8|. The QCD graphs that 
generate these distributions are called "Glasma graphs" . 
At high kx 3> Qs, the contribution of these graphs is 
highly suppressed. In contrast, for kx < Qs, where high 
occupancies in hadron wavefunctions are probed, Glasma 
graphs are enhanced by a^^, a factor of ~ 10^ for typi- 
cal values of the probed QCD fine structure constant as ■ 
In nuclear collisions at ultra-relativistic energies, these 



^ In hydrodynamical models, these long range rapidity correlations 
are a consequence of the choice of initial conditions, wherein the 
initial transverse spatial profile of the energy density distribution 
is assumed to be the same at all rapidities. Though not widely 
appreciated, this choice corresponds to an assumption of strong 



long range correlations in the dynamics of multiparticle produc- 
tion at short transverse spatial distances. Only azimuthal cor- 
relations arc dynamically generated by the hydrodynamic equa- 
tions. 



(nearly) boost invariant configurations arc argued to pro- 
vide the dominant mechanism for muhi-particle produc- 
tion, and factorization theorems (to leading logarithmic 
accuracy in x) have been derived J9l4ll|. 

Though the bulk of multi-particle production is nearly 
azimuthally isotropic, it is not exactly so. As first noted 
in [1^ , based on the formalism in [ll|, |l3[ , Glasma graphs 
produce contributions that are collimated at relative az- 
imuthal separations of A$ « and A$ « tt. It has 
been shown recently that these initial state contributions 
provide a quantitative description of the measured colli- 
mated yield in both proton-proton and proton-nucleus 
collisions 14l-ll7l|. 

However, as observed previously 18J, llSf, long range 
rapidity correlations from the initial state can also be 
collimated by the radial flow of a fluid. Indeed, within 
the Glasma flux tube framework itself, the radial flow of 
Glasma flux tubes correlated over distance scales 1/Qs 
was shown to generate a sizable ridge for large radial 
flow velocities [7|,|20|. In nucleus- nucleus collisions, where 
large radial flow is generated, several groups have shown 
that hydrodynamical flow provides a very good explana- 
tion of the data on two-particle correlations in the Ary- 
A$ plane |2ll424| . There have also been attempts to ex- 
tend this description of the ridge in nucleus-nucleus col- 
lisions to the ridges observed in high multiplicity proton- 
proton [H,!!!! and p-hPb ^^^ collisions. In the latter 
case, it is claimed that features of LHC high multiplicity 
data on proton-nucleus collisions [31[ and corresponding 
data in deuteron-gold collisions at RHIC [5| are quan- 
titatively explained in the Monte-Carlo (MC) Glauber 
hydrodynamic model of [27H29| . 

We will argue here that the applicability of hydrody- 
namics to the smaller size systems of proton-proton and 
proton/dcuteron-nucleus collisions is strongly dependent 
on assumptions about the nature of the initial multi- 
particle dynamics, much more so than in collisions of 
heavy nuclei. We will illustrate this by comparing results 
obtained in MC-Glauber models with particular dynam- 
ical assumptions about the initial state geometry with 
those obtained in the framework of the IP-Glasma ini- 
tial state model 



33| of hadrons and nuclei. Very 



noticeable differences are seen between the two models 
(with the same initial state configurations) for the com- 
puted eccentricities and corresponding flow coefflcients. 
In contrast, both initial state models, when combined 
with e vent-by-event hydrodynamical simulations, as in 
give similarly good descriptionfjj of bulk mul- 
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tiplicity and flow observables in heavy-ion collisions at 
both RHIC and the LHC. 



The paper is organized as follows. In the next sec- 
tion, we will outline the different methods employed to 
compute the initial spatial sizes and eccentricities and 
some of the consequences thereof. We will review the 
IP-Glasma model, and show its predictions for the initial 
spatial sizes in proton-proton and proton-nucleus colli- 
sions. We will compare the eccentricities obtained in this 
model to those in various implementations of the MC- 
Glauber model for proton-nucleus. The generated flow 
in proton-proton and proton/deuteron-nucleus collisions 
is considered next and contrasted between the two mod- 
els. The final section discusses the magnitude of viscous 
effects in different implementations of viscous hydrody- 
namics in proton-nucleus and nucleus-nucleus collisions. 
We end with a brief summary and outlook. 



MODELS OF THE INITIAL STATE GEOMETRY 

Modeling the initial state in p-l-A, d+A and especially 
p-|-p collisions is a lot more challenging than in A-f A col- 
lisions. In the latter, the system's geometry is primarily 
characterized by the overall shape of the interaction re- 
gion. The dominant component in shape fiuctuations are 
due to geometrical fiuctuations of nucleon positions in- 
side the nuclei from event-to-event. The large number of 
participants allows one, to first approximation, to neglect 
the dynamical details of how energy is deposited in A-l-A 
interactions. In p-l-A and d-|-A collisions, the system's 
geometry is very sensitive to the proton (or dcuteron) 
size, and the detailed nature of multi-particle production 
and the spatial distribution of the produced energy den- 
sity become important. In particular, sub-nucleon size 
fiuctuations (with characteristic length scales less than 
1 fm) contribute significantly to the initial geometry of 
matter produced in the collision. 




^ The IP-Glasma-|-MUSlC model of 1381 also reproduces the event- 
by-event v„ fluctuations measured by the ATLAS collabora- 
tion [33; at present, it appears to be the only model that suc- 
cessfully reproduces these flow fluctuations. 



FIG. 1. (Color online) Various models of the energy density 
deposition (denoted by red dots) in nucleon-nucleon collisions. 
In the left plot (a) the energy density is produced at the cen- 
ter of the colliding nucleons even for grazing collisions. The 
center and right plots (b),(c) correspond to different eccen- 
tricities depending on the matter distribution in the nucleon 
overlap region. For the configuration depicted on the left ec- 
centricity £2 ~ 1, whereas for the configuration in the center 
£2 =0. 



The spatial eccentricities that characterize the geome- 



try of the initial state can be defined as 



■\/(r" cos{ti(J)))'^ + (r" sm{n<j>))'^ 



where (•) is either an average over ah participant nucleon 
positions characterized by the nucleon centers or an av- 
erage weighted by the deposited energy density. 

In each realization of the MC-Glauber model, all par- 
ticipants contribute equally to the energy density de- 
posited in the system. For example, two nucleons that 
barely touched each other in the collision are assumed to 
deposit the same amount of energy as two nucleons that 
interacted with zero relative impact parameter. How- 
ever, in a microscopic parton model based picture of the 
hadron, it is natural to expect that peripheral nucleon- 
nucleon collisions deposit significantly less matter than 
the central ones. This is because the parton density 
decreases rapidly with impact parameter [40|, a picture 
confirmed by phenomenological descriptions of diffractive 
deeply inelastic scattering data from the HERA collider 

This 



in the framework of the IP-Sat dipole model [41 



physics is naturally incorporated in the IP-Glasma ap- 
proach, with the saturation scale depending on the im- 
pact parameter [32| . 

In many Glauber model computations, the energy even 
in a peripheral collision is deposited in the center of the 
wounded nucleons [42|. This is sketched in the left most 
configuration shown in Fig. [1] and corresponds to the 
ellipticity £2 ~ 1 and triangularity £3 = 0, when partici- 
pant centers are used to compute the average in Eq. ([T|) . 
However, on the basis of the parton model arguments 
outlined, one expects that peripheral nucleon-nucleon 
collisions will deposit most of the produced energy in 
the region of overlap, as illustrated in the middle figure 
of Fig. [TJ In this case, both the ellipticity £2 = and the 
triangularity £3 = 0, if one assumes an isotropic energy 
density deposition. Following the shape of the overlap 
region for the deposition of energy density, as sketched 
in the right most figure in Fig. [1] can produce instead 
eccentricities that are anywhere in between the two ex- 
tremes of the other two configurations. 

In A-l-A collisions, because of the large number of over- 
lapping nucleons, these finer details of geometry are less 
important. In contrast, these different microscopic pic- 
tures will have significant consequences for multiparticle 
production in nucleon-nucleon and nucleon-nucleus colli- 
sions. We note that in the IP-Sat framework, where par- 
ticle production only occurs in the geometrical overlap 
regions, a good description is obtained of the n-particle 
multiplicity distributions in proton-proton collisions at 
central rapidities from 200 GeV to 7 TeV [4i,|44|. 

We will now outline the IP-Glasma mo del, which has 
been discussed at length elsewhere [32, l33|. In this 



computed from fluctuating color charges via the classical 
QCD Yang-Mills equations. In an individual nucleon, 
the color charges are assumed to follow a local Gaussian 
distribution with variance (7^/Xp(a;,XT), a quantity pro- 
portional to the saturation scale Qs (^^i^t) h^ the pro- 
ton. This latter quantity is determined from the IP-Sat 
dipole model [4l| with parameters fit to HERA data on 
inclusive and exclusive final states [45| . Products of the 
nucleon dipole S-matrices, generate lumpy configurations 
of glue in nuclei, and are in agreement with extant fixed 
target data on electron- nucleus scattering [46| . This gives 
the variance .g^/i^ of Gaussian distributed nuclear charge 
distributions. The latter, in the MV model 43, IH], is 
used to solve for the coherent classical gauge fields of 
the nuclei before the collision. Here, as opposed to the 
previous implementation described in [32|, |33| , we com- 
pute X = Qs{x,XT)/y/s self-consistently at every trans- 
verse position x^. Further, the running coupling is evalu- 
ated locally as as{Qs{x,XT)), where Qs = nia.x{Qg,Qg) 
where Qg wee the saturation scales of the projectile 
(target). 

The solution for the transverse (longitudinal) gauge 
fields after the collision is obtained in Schwinger gauge 
A^ = and is given by the sum (commutator) of the 
incoming purely transverse fields to determine the initial 
energy and number content of the Glasma fields (4 
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The gauge fields are regulated by an infrared mass scale 
771 = 0.1 GeV, which ensures the unphysical Coulomb tail 
from solutions of the Yang-Mills equations is suppressed 
at large distances. This procedure, and the fact that the 
incoming fields are pure gauges, ensures that no energy 
density will be deposited in regions where at least one 
of the incoming fields vanishes. This key property of 
the solutions determines that in p+A (d+A) collisions 
the system size is dominated by the size of the incoming 
proton (and neutron). 

In Fig. [21 the system size rmax, the maximal radius for 
which the energy density of the Yang-Mills fields is above 
a minimal value of £ = q;Aq(-.q (with a G {1, 10}), is pre- 
sented. We observe that the system sizes in p-f-p and 
p-|-Pb are comparable and grow approximately linearly 
as a function of the number of gluons to the power of 
1/3. Rather than the difference between p-|-p and p+Pb, 
the result is more sensitive to the choice of the value of 
a, albeit relatively weakly given that the energy density 
changes by an order of magnitude. This uncertainty rep- 
resents intrinsic non-perturbative effects that cannot be 
further quantified within our present knowledge of QCD. 

Our observation of the scaling of the system sizes in 
p+p and p+Pb collisions indicates that their HBT radii|^ 



model, which goes si gnifi cantly beyond the treatment of 
hadron collisions in 43J, [4^, incoming gluon fields are 



^ It is interesting to note that a very similar trend to that of our 
^max is observed in p-|-p data on HBT radii measured at LHC 
[5511 . where all HBT radii grow approximately linearly with the 
number of produced particles to the power of 1/3. 




(dN/dy) ' 



FIG. 2. (Color online) System size in p+p and p+Pb collisions 
as a function of the number of gluons to the power of 1/3, 
{dN/dyY'^, computed in the IP-Glasma model, rmax is the 
maximal radius for which the energy density of the Yang-Mills 
fields are above a minimal value of Emin = aAgcjj, with a = 1 
(filled symbols) and a = 10 (open symbols). Note different 
energies for p+p and p+Pb. 



should be comparable in value. Indeed, for the same 
number of produced particles and comparable sizes of two 
systems, leading to comparable energy densities, both 
systems are very similar and we do not expect to observe 
significant differences in their HBT radii. Let us note 
that we do not calculate explicitly HBT radii in p+p and 
p+Pb, we only compare both systems based on their sizes 
and energy densities. If subsequently fluid dynamical 
evolution in p+A collisions is significant relative to p+p 
collisions, we would then anticipate significantly different 
HBT radii in p+Pb, as recently discussed in Ref. J29| . 
Thus HBT radii will help to discriminate between models 
of the spatial distribution of matter in the initial state 
and the magnitude of radial flow experienced in each. 

In Fig. m the ellipticity and triangularity computed in 
the IP-Glasma model are plotted as a function of dN/dy, 
the gluon number per unit rapidity in the model. The 
ellipticities in p+p collisions are significantly lower than 
in p+Pb collisions, except at very low values of the mul- 
tiplicity. The triangularity £3 in p+p collisions is very 
small and is comparable to the ellipticity. The triangular- 
ity in p+Pb collisions is consistently larger than in p+p, 
though distinctly smaller than the ellipticity in p+Pb up 
to very high multiplicities. 

In Fig. m we plot £2, £3 in proton- lead collisions in 
the IP-Glasma model as a function of the number of par- 
ticipants A^part and compare the results to two different 
realizations of the Monte-Carlo Glauber model. Com- 
puting eccentricities using participant centers in a MC- 
Glauber model, as done in |27| , the results are along the 
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lines anticipated in our discussion of Fig. [T] £2 is exactly 
unity and £3=0 for A^part — 2. A similar trend for £2 at 



FIG. 3. (Color online) Eccentricities £2 and £3 as a function 
of dN/dy in the IP-Glasma model for p+Pb (filled symbols) 
and p+p (open symbols) collisions. 



low A^part is observed when computing eccentricities us- 
ing Gaussian energy densities (with width ctq = 0.4 fm) in 
the centers of participants (MC-Glauber 1) as a weight. 
In this case the eccentricities are noticeably lower. The 
MC-Glauber 1 realization of the Glauber model is similar 
to the one employed in the computations of 27H29| for 
proton-nucleus and dcuteron-nucleus collisions. 

Assigning a Gaussian distributed energy density to the 
midpoint between two colliding nucleons (MC-Glaubcr 
2), as illustrated in the middle figure in Fig.[Tl reduces the 
resulting eccentricities significantly, with the difference in 
£2 between models 1 and 2 being approximately a factor 
of 2. We have checked that decreasing the smearing width 
from CTo = 0.4 fm increases both eccentricities. 

Wc have used a black disk approximation of the cross- 
section, meaning that nucleons arc wounded whenever 
their geometric distance from a nuclcon of the other nu- 
cleus is less than r^vAr = \/<Jnn/tt, where aNN is the 
nuclcon-nuclcon inelastic cross section. Alternatively, 
one can introduce a smooth profile of the nucleon that 
determines the probability for an interaction at a given 
nucleon-nucleon distance. This profile can be extracted 
from the p+p differential elastic cross section data [56| 
and can be approximated by a Gaussian 57H59l|. Its use 
has been argued to be preferable because it does not re- 
sult in extremely large elastic nuclcon-nuclcon cross sec- 
tions as the hard-sphere case. Using a smooth profile, we 
find an increase in both the system size (by up to 50%) 
and the eccentricities (by up to a factor of 2). Again, 
it demonstrates that eccentricities in p+Pb collisions are 
very sensitive to details of nucleon-nucleon interactions. 

In the IP-Glasma model, both the ellipticity and tri- 
angularity coincide with model MC-Glauber 2 for all but 
the smallest values of A^part- This agreement is however 
a coincidence for the value of (Tq = 0.4 fm chosen in the 
MC-Glaubcr model and will not hold if this parameter 
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FIG. 4. (Color online) Eccentricities e2 (a) and £3 (b) as 
a function of the number of wounded nucleons A'^part . In the 
MC-Glauber (participant centers) model the energy density is 
deposited in the centers of wounded nucleons (without smear- 
ing). Smearing energy densities with the Gaussian distribu- 
tion ((To = 0.4 fm) results in the MC-Glauber 1 model. In the 
MC-Glauber 2 model the energy density is smeared about the 
midpoint between colliding nucleons. 



obtain the thickness function Tdirr)- 

In Fig. [5l we show plots with a typical deuteron config- 
uration in the transverse plane (denoted by open circles) 
superposed on transverse projections of the nucleon posi- 
tions in the gold nucleus (denoted by filled circles). In the 
top plot, we show the energy density contours from the 
MC-Glauber 1 model and in the bottom plot, the corre- 
sponding IP-Glasma model results. These are seen to be 
quite different. In the latter, it is observed that the peaks 
in the contour are closely associated with the centers of 
the deuteron nucleon positions and vary strongly depend- 
ing on the number of gold nucleon positions in their im- 
mediate vicinity. In the former MC-Glauber case, signif- 
icant energy densities are seen even in regions where nu- 
cleons of the gold nucleus are widely separated in trans- 
verse spatial position from the deuteron nucleons. Nucle- 
ons that have been marginally grazed produce as much 
energy density as those that have suffered a head on colli- 
sion. In the IP-Glasma model, because the mean distance 
in the projection onto the transverse plane between the 
two nucleons in a deuteron is 2.52 fm, the majority of 
events have widely separated interaction regions. This is 
quite different in the MC-Glauber model. 

Whether eccentricity is a relevant measure in deuteron- 
gold collisions depends sensitively on the radial separa- 
tion of the regions where energy density is deposited. If 
they are too far apart for hydrodynamic flow to bring 
them into contact over the system's lifetime, the eccen- 
tricity will be a poor measure of flow. If they are close 
enough at the same eccentricity to influence subsequent 
flow, the eccentricity will track flow better. Thus eccen- 
tricity in deutcron-gold collisions, in contrast to nucleus- 
nucleus collisions, is at best a qualitative measure of 
anisotropic flow. 



FLOW IN PROTON-PROTON AND 
PROTON/DEUTERON-NUCLEUS COLLISIONS 



is varied. In the IP-Glasma model, there is no such free 
parameter. 

The differences between the MC-Glauber realizations 
used in hydrodynamical models and the IP-Glasma 
model are strikingly seen in deuteron-gold collisions. For 
the deuteron, we use the Hulthen form of the wave func- 
tion. 



4>pn{r) 



1 ^J ah{a + h) e 



-hr 



(2) 



with r being the distance between the proton and the 
neutron and the parameters a = 0.228 fm"^ and b = 
1.18 fm~^ [60|. Integrating cj)^^ over z = ^,.2 _ ^2^^ 
where r^ is the distance between the proton and the 
neutron in the projection onto the transverse plane, we 



An interesting compilation of the ratio of the elliptic 
flow coefficient to the ellipticity V2/S2 versus the mul- 
tiplicity at central rapidity in proton- nucleus, deuteron- 
nucleus, and a variety of centralities in nucleus-nucleus 
collisions can be found in |5| . This scaling may be taken 
to suggest the same pattern of collective flow in the 
smaller size systems as in the larger size one. However, 
as noted, the effect of the different eccentricity compu- 
tations on vile2 for p/d-A collisions is dramatic because 
£2 in the IP-Glasma model is two to five times smaller 
than the MC-Glauber model depending on A^part- For 
A+A collisions, the differences in the eccentricity com- 
putations are not as large [3^, and the scaling of W2/£2 
is less sensitive to model assumptions. 

In addition, from the point of view of examining the 
presence of collective flow in the system, a more useful 
variable to plot v^js^ against is the multiplicity per unit 




FIG. 5. (Color online) Initial energy density distribution (ar- 
bitrary units, increasing from blue to red) in the transverse 
plane in a d+Au collision in the MC-Glauber model (a) and 
the IP-Glasma approach (b). The nucleon positions (open 
circles for the deuteron, solid circles for gold) are exactly the 
same in the two cases. 



transverse overlap area S± (61[. Such a compilation, for 
a wide range of centralities in heavy-ion collisions at the 
LHC and RHIC can be found in [62| . With this criterion 
however, when the x-axis is divided by S±, the scaling of 
V2/S2 is broken immediately for d+Au collisions. Indeed, 
the initial transverse area in d+Au is approximately a 
factor of two larger than in p-f-A collisions. It will be in- 
teresting to calculate the transverse area for all colliding 
systems systematically in different models and investi- 
gate further the model dependence of scaling of V2/S2 in 
smaller size systems. This study will be reported else- 
where. 

To study whether the behavior of the eccentricities pre- 
sented in Fig. U is a good representation of the generated 



flow if the system behaves hydrodynamically, we com- 
pute the root mean square V2 and V3 integrated over 
Pt > 0.5 GeV for different TVpart in the MC-Glauber 1 
and IP-Glasma model. We use a constant 77/s = 0.08, an 
initial time of hydrodynamic evolution tq = 0.2 fm and 
a freeze-out temperature of Tfo = 120 MeV. The nor- 
malization of the initial energy density in the IP-Glasma 
model was tuned to reproduce the charged particle mul- 
tiplicity in p-l-p collisions at 7 TeV. In the MC-Glauber 1 
model the normalization of the energy density was set to 
approximately produce the same amount of charged par- 
ticles as the IP-Glasma model. Because wc are only in- 
terested in general trends in this work and not a detailed 
comparison to experimental data, we have not performed 
any fine tuning of parameters to reproduce particle spec- 
tra in p/d-|-A collisions. The reader should note however 
that these initial conditions, specifically the very low rj/s 
and small tq, can reasonably be considered to provide 
upper bounds on the magnitude of the generated flowO 

We show results for the integrated root mean square 
V2 and V3 for p-|-Pb collisions in Fig. [6] and for d-|-Au 
collisions in Fig. [71 The first thing to note is the qualita- 
tive difference between the ccntrality dependence of V2 in 
p+Pb and d+Au collisions. While in p+Pb collisions V2 
drops with increasing iVpart as expected from e2(-^part) 
in d+Au we find the opposite behavior. This behavior is 
expected qualitatively from e2(^part) in d+Au collisions 
27|. However, as per our discussion in the previous sec- 
tion, £2 alone is not necessarily useful for a quantitative 
understanding of flow in d+Au collisions. 

In p+Pb collisions, given the eccentricity £2 of the 
MC-Glauber 1 model in Fig.[H one might naively expect 
an increase of V2 by a factor of three when going from 
^part = 20 to A^part = 7 if it scalcs with £2- While we do 
find the same trend, V2 changes by a relatively smaller 
factor of approximately 1.7. In the IP-Glasma model, the 
change in V2 with A^part is larger (a factor of 2.5) even 
though the eccentricity £2 varies more slowly than in the 
MC-Glauber 1 model. For Apart = 14, V2 is approxi- 
mately a factor of two smaller than in the MC-Glauber 1 
model, and about 60% lower for Apart = 20. ^3 is nearly 
flat in the MC-Glauber 1 model and decreases with Apart 
in the IP-Glasma model for both p+Pb and d+Au colli- 
sions. 

We conclude that for small size systems, like p+Pb 
or d+Au, there is no simple quantitative scaling of the 
flow with eccentricity. Further, a smaller ellipticity in a 
different initial state model does not necessarily lead to 
smaller V2 in that model, because the geometries (and 
system sizes) may be so drastically different, that £2 is 
not a sufficient predictor of V2. This is seen strikingly 



* For nucleus-nucleus collisions at RHIC and the LHC, average 
values of rj/s = 0.12 and r]/s = 0.2 respectively, give the best fits 
to data [631. 



for the flow generated in p+p collisions relative to p+Pb 
collisions. On the basis of the plots in Fig. |3l one might 
conclude that flow is much smaller in p+p relative to 
p+A. 

We computed the integrated {px > 0.5 GeV) 
anisotropic flow for p+p collisions at y^ = 7 TeV in the 
IP-Glasma model. We find {v^y/^ « 0.02 at b = Ofm 
and (w|)^/2 ~ 0.035 at & = 1 fm - for 6 = fm the muhi- 
plicities in p+p are typically large, and the results can be 
qualitatively compared to those in p+Pb collisions. As 
shown in Fig. [HI W2 at A'^part = 20 in p+Pb is comparable 
to the results in p+p collisions within 50%. The lat- 
ter values, computed in the IP-Glasma model, arc given 
by the points at A^part = 2 in Fig.[6l We further find 
(f|)^/^ ~ 0-01 for both studied impact parameters in 
p+p collisions, to be compared to the values shown for 
p+Pb in Fig.El 

We next present V2 and v^ as functions of transverse 
momentum px for p+Pb collisions in Figs. [8] and |9] and 
for d+Au collisions in Figs.lTOlandfTTl We see the same 
trend as observed for the integrated t;„. 
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FIG. 6. (Color online) Integrated {v1)^^^ and (^1)^'^^ for 
charged hadrons in p+Pb collisions at different A'part in the 
MC-Glauber 1 and IP-Glasma model for pr > 0.5 GeV and 
ri/s — 0.08. V2 decreases with A'^part. Results for p+p colli- 
sions are for 6 = fm. 



Finally, we can qualitatively compare our results here 
to the results of the LHC and RHIC experiments on 
proton-nucleus and dcuteron-nucleus collisions. The 
trend of V2 as a function of ccntrality observed in p+Pb 
collisions appears to be different from that of the AL- 
ICE data [3| on proton-lead collisions at y^ = 5.02 
TeV/nucleon. However, the error bars in the pub- 
lished data are too large to draw a definitive conclu- 
sion at present. The ATLAS collaboration has also pre- 
sented [3l[ a quantity called 112 (PC), which is defined 
similarly to the ALICE V2 and has the same trend as the 
ALICE results. However, the collaboration also presents 
results for t!2{4}, from four particle correlations, which 
appears to have the opposite trend with centrality rela- 
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FIG. 7. (Color online) Integrated {v^)'^^^ and (wf)^''^ for 
charged hadrons in d+Au collisions at different A''part in the 
MC-Glauber 1 and IP-Glasma model for pr > 0.5 GeV and 
77/5 — 0.08. V2 increases with A''part. 
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FIG. 8. (Color online) V2{pt) for charged hadrons in p+Pb 
collisions at fixed A'part = 7 (thin lines) and 20 (thick lines) 
in the MC-Glauber (dashed) and IP-Glasma (solid) model. 



five to W2(PC). The computations of [30| appear to be 
in agreement with this W2{4} quantity for the centralities 
compared. However, as can be seen in Fig. [HI the IP- 
Glasma results are approximately a factor of two lower 
for the Apart that correspond to the same centrality se- 
lection. Note further that the IP-Glasma results are for 
ij/s = 0.08 and will be smaller for the rj/s = 0.2 that 
gives a good description of V2 in A+A collisions at the 
LHC. We also note that the RHIC d+Au results on the 
ridge are reproduced by the MC Glauber 1 model. The 
differences between this model and the IP-Glasma model 
(for integrated V2 values) in these collisions can be seen 
in Fig. [71 More quantitative studies and additional data 
will clearly help clarify the role of hydrodynamics in the 
interpretation of the RHIC and LHC results on the ridge 
in dcuteron-gold and proton-lead collisions respectively. 
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FIG. 9. (Color online) W3(pt) for charged hadrons in p+Pb 
collisions at fixed A'^part ~ 7 (thin lines) and 20 (thick lines) 
in the MC-Glauber (dashed) and IP-Glasma (solid) model. 



FIG. 11. (Color online) D3(pt) for charged hadrons in d+Au 
collisions at fixed A'^part = 10 (thin lines) and 30 (thick lines) 
in the MC-Glauber (dashed) and IP-Glasma (solid) model. 
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FIG. 10. (Color online) 112 (pr) for charged hadrons in d-|-Au 
collisions at fixed TVpart = 10 (thin lines) and 30 (thick lines) 
in the MC-Glauber (dashed) and IP-Glasma (solid) model. 



VALIDITY OF VISCOUS HYDRODYNAMICS 

In the previous sections, we discussed the strong de- 
pendence of the initial geometry on initial state dynamics 
when one considers the especially small sized systems one 
has in proton/deutcron-nuclcus collisions. We also dis- 
cussed the hydrodynamic flow resulting from these initial 
spatial geometries. One may ask in addition under what 
conditions, if any, hydrodynamics is applicable to these 
especially small sized systems. Hydrodynamics is usually 
a good effective field theory in the late time, long wave- 
length limit of the theory. One way to quantify whether 
this holds in the systems of interest is to compare the rel- 
ative magnitude of the viscous terms in the stress-energy 
tensor to the ideal fluid terms. Considering for simplic- 
ity only shear effects, and neglecting heat flow and bulk 



rji^f rpllU 



(3) 



where the ideal term Tl^" == (e + P)g^"' - Pg^''' and the 
viscous part of the stress-energy tensor satisfies the equa- 
tion 



t,A/A/i.^9a^"^ + TT^'' = --^^•'{dxu'') 



V 






(4) 



Here, e is the energy density and P the pressure, the 
metric g'^''' = diag(l, —1, —1, —1), the transverse projec- 



tor A^" 



9 



fiiy 



uf^u" and V^ ^ At'^d^. The relaxation 



time is set to r^ = 877/(6 -|- P). 

We cannot initialize with the full T^" provided by the 
IP-Glasma model at the initial time, because of its highly 
non-equilibrium nature. Therefore, we are left with a 
choice for ttq'^ at the initial switching time r = 0.2 fm 
between flow in the Glasma and the later flow described 
by hydrodynamics. We study the case ttq" = 0, as often 
employed in viscous hydrodynamic simulations, and the 
case where ttq'' takes on its Navier-Stokes value. 



^ V 



Vu" - -A'^'Vxu^ 



(5) 



In this work wc neglect gradients of the transverse veloc- 
ities at the initial time and only keep the dominant piece 
from the longitudinal (boost-invariant) dynamics. 

To avoid instabilities of the algorithm outside the 
freeze-out surface, where viscous corrections can become 
very large, we introduce a regulator that restricts the 
viscous correction to be smaller than ten times the ideal 

< 



part. In practice, we require y^Tr^I^/TfP^ < lOVe^ + 3P^ 
and implement a continuous regulation of tt'^'^ to satisfy 



this requirement: 



^f 



:.liv 



3P2). 



= ■K^'^ tanh(p)/p , where 
We have checked that the 



results presented for ccUs within the freeze-out surface 
are very weakly sensitive to this regulator when the code 
is stable without regulation. The only effect of the reg- 
ulator is to avoid instabilities that are triggered outside 
the freeze-out surface. 

To quantify the validity of the viscous hydrodynamic 
approach for different collision systems, wc can deter- 
mine the fraction of cells within the freeze-out surface 



yTT^"^ that are either 



that have a viscous correction y^j,y/ 
larger than 25% (Figs. [T^ and [T^ of the ideal fluid con- 
"3P2^ or larger t han a 50% (Fig. ^ to 
3P2 plays the role 



tribution ve^ 

the same. The ratio ^/ttJIuT^ /\e^ 

of an effective inverse Reynolds number (6^, |6£ 



The results shown are for two different values of rj/ s, 
0.08 and 0.2 and arc for central (6 = Ofm) Pb+Pb and 
central p-f Pb collisions as a function of time r for zero 
initial tt^'' in Fig. [121 and for an initial Navier-Stokes 
TT^'' in Figs. [13] and [T4| The results are averages over 
10 events each. Wc find very similar behavior in p+Pb 
and Pb-|-Pb collisions, with many cells having large cor- 
rections at times t < 2 fm. especially for the larger value 
of rj/ s — 0.2. In the latter case, which gives the best 
agreement with heavy ion collision data at the LHC [63| , 
nearly all cells have corrections of at least 25% at early 
times, even for zero initial tt^''. In the case of Navier- 
Stokes initial conditions, all cells within the freeze-out 
surface start out with an at least 50% viscous correction. 
For the smaller 77/s = 0.08, viscous corrections are small 
after t = 1 fm/c for both initial ir^" choices. 

Although the results seem very similar for p+Pb and 
Pb-f Pb collisions, it is important to note that the lifetime 
in Pb-I-Pb collisions is about 6 times longer than in p-|-Pb 
collisions. This means that viscous corrections are large 
for a significant fraction of the total lifetime of p-|-Pb col- 
lisions, while they are large only for a small fraction of 
the total lifetime of Pb-I-Pb collisions. In other words, 
depending on the value of 77/s, viscous hydrodynamics in 
the first 1-2 fm/c is sensitive to not just the initial spa- 
tial geometry but also to details of how the viscous flow 
is initialized. While there is still a significant part of the 
space-time evolution where viscous hydrodynamics is a 
reliable framework for Pb-|~Pb collisions, the same can- 
not be said for p-|-Pb collisions. In this case, the system 
dilutes by 3-4 fm/c - thus, viscous hydrodynamics is un- 
reliable for approximately half the lifetime of the system 
in p-|-Pb collisions, and therefore a source of potentially 
large systematic errors. Similar conclusions and limits 
on the initial time tq when viscous h ydr odynamics be- 
comes applicable have been found in [65| for heavy-ion 
collisions. 
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FIG. 12. (Color online) Fraction of all cells within the freeze- 
out surface that have a viscous correction of at least 25% of 
the ideal T^"" . Initialization with tt^" = 0. 
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FIG. 13. (Color online) Fraction of all cells within the freeze- 
out surface that have a viscous correction of at least 25% of 
the ideal y". Initialization with the Navier-Stokes value for 



SUMMARY AND OUTLOOK 

In this work, we have demonstrated that different 
dynamical assumptions can lead to qualitatively differ- 
ent results for spatial ellipticities and triangularities in 
smaller size systems. In contrast, the same range of 
dynamical assumptions give qualitatively comparable re- 
sults in heavy ion collisions, with some differences show- 
ing up in quantitative studies. Results for the flow co- 
efficients V2 are significantly different between the MC- 
Glauber 1 model and the IP-Glasma model in proton- 
nucleus collisions at the LHC. The results studied are 
for 77/s = 0.08, and one expects that V2 in both models 
will be significantly smaller for the value of 77/s = 0.2; 
this value best describes the LHC Pb-fPb data in the 
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FIG. 14. (Color online) Fraction of all cells within the freeze- 
out surface that have a viscous correction of at least 50% of 
the ideal T^" . Initialization with the Navier-Stokes value for 



IP-Glasma model. The trend of V2 with centrality seen 
in both models is different from the preliminary ALICE 
V2 and ATLAS U2(PC) data on p+Pb colhsions. More 
precise data from all the LHC experiments can resolve 
this issue~the four particle correlations measured by AT- 
LAS are an interesting step in this direction. With re- 
gard to deuteron-gold collisions at RHIC, the differences 
in V2 between the MC-Glauber 1 and IP-Glasma model 
are smaller; again, both results would be smaller for the 
Tj/s = 0.12 value that gives the best fit to RHIC Au-|-Au 
data in the IP-Glasma model, v^^ even for the small rj/s 
used here, is quite small at both RHIC and the LHC. 

In addition to details of the initial spatial geometry, 
flow in small size systems can have large viscous contri- 
butions for a significant fraction of the lifetime of the 
system. Within the framework of the IP-Glasma-|-MUSiC 
model, we showed for instance that for rj/s = 0.2, viscous 
corrections are as large as 50% for times up to a fermi. 
In the Pb-|-Pb case, much of the flow is built up at later 
times when these viscous corrections are small, and sec- 
ond order viscous hydrodynamics can be considered a 
good effective description of the collective dynamics. In 
contrast, because the lifetime of the system in p-|-Pb col- 
lisions is much shorter than in Pb+Pb, the large viscous 
corrections over a signiflcant fraction of the lifetime sug- 
gest second order viscous hydrodynamics is less reliable 
in such systems. 

Hydrodynamics has provided us with a very success- 
ful framework to interpret the results of momentum 
anisotropies observed in heavy ion collisions. The recent 
experiments on proton/deuteron-nucleus collisions have 
brought to the fore the question of the applicability of 
hydrodynamics to these very small sized systems. The 
interest in further studies along the lines of this work 
is two-fold. Firstly, it would be important to quantify 



where hydrodynamics breaks down to establish a better 
understanding of the transport properties of the Quark- 
Gluon plasma. Here a careful study of the ratio U2/e2 as 
a function of the transverse density, Nch/S±, in various 
systems could be very useful. However, as we emphasized 
in this paper this problem is not trivial because in small 
systems various models lead to very different eccentrici- 
ties. 

Secondly, these studies focus our attention on a better 
understanding of the initial non-equilibrium dynamics of 
strong color fields. For instance, if the azimuthal collima- 
tion of the nearside p+p, p+Pb and d+Au ridge is not 
due to hydrodynamic flow, a subtle long range quantum 
interference between saturated gluon fields is likely the 
successful alternative explanation. On the other hand, if 
hydrodynamics is a viable explanation, quantitative com- 
parisons to data can help narrow down the dynamics that 
generates the initial spatial geometries that subsequently 
experience hydrodynamic flow. 
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